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Abstract. - In this work we discuss if fermionic sources could be responsible for accelerated 
periods in a Friedmann-Robertson- Walker spatially flat universe, including a usual self-interaction 
potential of the Nambu-Jona-Lasinio type together with a fermion-scalar interaction potential of 
the Yukawa type. The results show that the combination of these potentials could promote an 
initially accelerated period, going through a middle decelerated era, with a final eternal accelerated 
period, where the self-interaction contribution dominates. 
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The identification of constituents that promoted the 
inflationary period in the evolution of the universe is 
a fundamental topic in cosmology. Several candidates 
has been tested for describing both the inflationary pe- 
riod and the present accelerated era: scalar fields, exotic 
equations of state and cosmological constant [T] . Besides 
those, fermionic fields has also been tested as gravita- 
tional sources of an expanding universe. In fact, fermionic 
sources can be responsible of accelerated periods with dif- 
ferent regimes emerging from it [2H5]- In some of these 
models the fermionic field plays the role of the inflaton 
in the early period of the universe and of dark energy for 
the old universe, without the need of a cosmological con- 
stant term or a scalar field. In an old universe scenario an 
initially matter dominated period gradually turns into a 
dark (fermion) energy period when an accelerated regime 
starts and remains for the rest of the evolution of the sys- 
tem. These fermionic sources has been investigated using 
several approaches, with results including numerical so- 
lutions, exact solutions, anisotropy-to-isotropy scenarios 
and cyclic cosmologies (see, for example [2H5])- When con- 
sidering these models, one important point is the choice 
of the fermionic potential and in previous works [3] self- 
interacting potentials were tested to account for acceler- 
ated regimes. One complementary/alternative approach 
would be consider interactions between a scalar and a 
fermionic field through the presence of a potential of the 



(a) gravitam@yahoo.com 

devecchi@fisica.ufpr.br 
f"^' kremer@fisica.ufpr . br 



Yukawa type [5]. This potential was proposed originally 
in particle physics to describe the behavior of strong force 
interactive fermions and here we want to test which kind 
of role the Yukawa potential could play in an young ac- 
celerating universe. Besides, as in previous works [2113], a 
fermionic self-interaction term of the Nambu-Jona-Lasinio 
type [B] is included. We have used the metric signature 
— , — , — ) and natural units with SttG = c = h = 1. 

Let us start with a brief review of the elements of the 
tetrad formalism employed to include fermionic fields in a 
dynamical curved space-time, since, as it is well known, 
the tetrad formalism permits the inclusion of fermions 
in gravitational models. Following the general covariance 
principle, a connection between the tetrad and the met- 
ric tensor g^i, is established through the relation = 
e°e^77afc, a = 0, 1, 2, 3 where ej^ denotes the tetrad or vier- 
bein and r]ab is the Minkowski metric tensor (see e.g. [ZHS]). 
Here Latin indices refer to the local inertial frame whereas 
Greek indices to the general system. Furthermore, the 
general covariance principle imposes that the usual Dirac- 
Pauli matrices 7° must be replaced by their generalized 
counterparts F^ = ej^7°, where these matrices satisfy the 
extended Clifford algebra, i.e., {F'',r''} = 2g'"'. 

The generally covariant Dirac Lagrangian density for 
the fermionic field is given by 

Cmr.ci^') = '^H'r^D^i: - {D^i^)T'^ij] - F(7^), (1) 

where and 1/' = ■ip^j'^ denote the spinor field and its ad- 
joint, respectively, and V^ip) is the self-interaction poten- 
tial of the fermionic field which is supposed to be massless. 
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In the above equation the covariant derivatives read 

Df,i^ = - n^i^, D^^ = d^^ + i^n^, (2) 

where the spin connection is given by 
1 



(3) 



with rj^_^ denoting the Christoffel symbols. 

The Lagrangian density for a massive scalar field 4> and 
the one corresponding to the Yukawa interaction between 
the fermionic and the scalar field read 



'CYukawa((/', "0) 
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where m is the mass of the scalar field and A the coupling 
constant of the Yukawa potential. 

The total action for a massless fermionic field and a mas- 
sive scalar field that are connected through an interaction 
of Yukawa type can be written as 



S 



1 



1 



1 



(6) 



where R is the scalar curvature. 

We obtain the Dirac equations for the spinor field and its 
adjoint coupled to the gravitational field from the Euler- 
Lagrange equations for ip and tp, namely 

dV_ 

dV - 



= 0, 



(7) 
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Likewise the Euler-Lagrange equation for </) leads to the 
modified Klein-Gordon equation 



(9) 



From the variation of the total action ^ with respect 
to the tetrad we obtain Einstein field equations 
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where T^i/ is the total energy-momentum tensor which is 
a sum of the contributions from the fermionic and scalar 
fields. Since we arc dealing with a fermionic field in a 
space-time without torsion, the total energy-momentum 
tensor is symmetric and reads 
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-AV'cAV + 2 (V'r^-DAV' - i^A^Z-rV) - v{ii;)\ . (11) 



The Friedmann-Robertson- Walker metric 

ds^ ^ dt^ - a{tf{dx^ + dy^ +dz^), (12) 

with a{t) denoting the cosmic scale factor, mirrors the 
homogeneity and isotropy properties of a spatially flat 
universe. For this metric the components of the tetrad, 
Dirac-Pauli matrices and spin connection read 



(13) 



a{t) 



where the dot represents the derivative with respect to 
time. 

Furthermore, from the hypothesis of homogeneity and 
isotropy it follows that the spinor and scalar fields depend 
only on time so that the Dirac (O, (|8]) and Klein-Gordon 
(HJ equations reduce to 



Z a dtp 

^ + ^ _ _ jA#7° = 0, 

2 a dtp 

(j) + 3-0 -I- rri^cf) + Xlpilj = 0. 
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From the Einstein field equations (|10p it follows the 
Friedmann and acceleration equations 
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(p + 3p). 



(18) 



respectively, where the total energy density p and total 
pressure p are are given by 

P=l^^ + l^^^^ + V(V^'A) + Vii;), (19) 

We suppose now that the self-interaction potential of 
the fermionic field is given by y = A i^ipi') , where A is 
a coupling constant and n a constant exponent. If we 
substitute this potential into the Dirac equations (|15p and 
(|T6| . multiply the first one by ip, the second one by ip 
and add the resulting equations we obtain a difl'erential 
equation for the bilinear ipip whose solution is given by 
tpip — C /a^ where C denotes an integration constant. 

Once we know the explicit form of the bilinear = 
Cjc? the total energy density (|T9| and total pressure (|20| 
become 
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where we have introduce the constants C\ = C\ and Ci ~ 
CA. 
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Fig. 1: Acceleration field a as a function of time t. 
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Fig. 2: Yukawa potential il)(j>ili as a function of time t. 



Therefore, from p7)) and ([T5)) 9 we have a coupled system 
of non-linear differential equations for the determination 
of a{t) and (j}{t) that include the modified Klein-Gordon 
equation and the acceleration equation, namely. 
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20^ - m^cf^ + 
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The coupled system of differential equations and 
p4|) was solved numerically with the initial conditions that 
the cosmic scale factor and the energy density were nor- 
malized at t = 0, i.e., a(0) = 1 and p{Q) = 1, which from 
the Friedmann equation ((T5)) i implies that d(0) = l/-\/3. 
Furthermore, it was supposed that the scalar field at time 
t = was zero with a very small slope, i.e., 0(0) = 
and (^(0) = 10~^. With these initial conditions, the cou- 
pling constant related with the self-interaction potential 
of the fermionic field is no more arbitrary, since from 
the expression for the energy density (PT|) we must have 
C2 = p{0) — 0(0)^/2. However, it remains to specify the 
coupling constant of the Yukawa potential Ci , the mass of 
the scalar field m and the exponent of the self-interaction 
potential of the fermionic field n. We have chosen fixed 
normalized values for the mass of the scalar field and for 
the exponent of the self-interaction potential of the scalar 
field - namely, m — 10"^ and n ~ 0.5 - and variable values 
for the coupling constant of the Yukawa potential in order 
to see the influence of the strength of this potential on the 
solutions of the coupled system of differential equations. 

As a result of the numerical integration, the scale factor 
behaves accordingly to an ever-expansive universe. The 
classification of eras can be done in terms of the time evo- 
lution of the acceleration field. In fact, in Figure 1 we 
have plotted the behavior of the acceleration a as a func- 
tion of time, for three different (normalized) values of the 
coupling constant of the Yukawa potential: |Ci | = 2.0, 2.4 
and 2.8. From this figure it is possible to infer that we have 



three periods which can be interpreted as an initial infla- 
tionary period where the acceleration is positive, a period 
dominated by radiation/matter where the acceleration is 
negative and a period dominated by the dark energy where 
the acceleration becomes again positive. Furthermore, we 
observe that longer periods of deceleration are related with 
larger values of |Ci|. By increasing the value of |Ci| be- 
yond 2.8 the solution of acceleration field presents only an 
accelerated period followed by a decelerated one, while for 
values of \Ci\ smaller than 2.0 the acceleration field has 
only one solution corresponding to an accelerated period. 

In Figure 2 we show the behavior of the Yukawa poten- 
tial i/jc/)^ as function of time for different values of the cou- 
pling constant Ci . The behavior of the curves in this figure 
corroborates with the explanation given above, i.e., by in- 
crease of the coupling constant the amplitude of Yukawa 
potential becomes more accentuated, so that it has a di- 
rect influence on the decelerated period. 

The coupling constant of the Yukawa potential has also 
an influence in the time decay of the total energy density. 
This fact can be observed from the analysis of Figure 3 
where the total energy density is plotted as a function of 
time for different values of Ci. 

To sum up, in this work we have investigated the pos- 
sibility of a fermionic cosmology to describe the different 
accelerated regimes of our universe. This is possible due 
to the inclusion of two potentials, one self-interactive and 
one of the Yukawa type. When the universe is still young, 
the combination of these contributions promote a short 
positive accelerated regime that can be associated to a in- 
flationary period, followed by a decelerated period that 
would correspond to a matter/radiation dominated uni- 
verse. We conclude that the intensity of the Yukawa cou- 
pling determines the time interval and the amplitude of 
the decelerated regime. On the other hand, the constant 
that controls the intensity of the self-interactive potential 
is decisive to the existence of a final accelerated period, 
that could correspond to a universe dominated by dark 
energy. 
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Fig. 3: Total energy density p as a function of time t. 
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